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Abstract 

In this paper we consider a mathematical model for the evolution and collapse of 
the Easter Island society, starting from the fifth century until the last period of the 
society collapse (fifteen century). Based on historical reports, the available primary 
sources consisted almost exclusively on the trees. We describe the inhabitants and 
the resources as an isolated system and both considered as dynamic variables. A 
mathematical analysis about why the structure of the Easter Island community 
collapse is performed. In particular, we analyze the critical values of the fundamental 
parameters driving the interaction humans-environment and consequently leading 
to the collapse. The technological parameter, quantifying the exploitation of the 
resources, is calculated and applied to the case of other extinguished civilization 
(Copan Maya) confirming, with a sufficiently precise estimation, the consistency of 
the adopted model. 
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1 Introduction 



Easter Island history is a very famous example of an evolved human society 
that collapsed for over exploiting its fundamental resources [1 .2.3J that in this 
case were essentially in palm trees. It were covering the island jl] when, few 
dozens of individuals, first landed around 400 A.D. Its advanced culture was 
developed in a period of one thousand years approximately. Its ceremonial 
rituals and associated construction were demanding more and more natural 
resources especially palm trees. The over exploitation of this kind of tree, 
very necessary as a primary resource (tools construction, cooking, erosion 
barrier, etc.) was related with the collapse. In this paper, a mathematical 
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model concerning growing and collapse of this society is presented. Different 
to usual like Lotka-Volterra models [5|6] where the carrying capacity variation 
becomes from external natural forces, in this work it is directly connected 
with the population dynamics. Namely, population and carrying capacity are 
interacting dynamics variables, so generalizing the Leslie model prey-predator. 

The general mathematical treatment of a model describing such a complex 
society is a very hard task and probably not unique. Our aim is to settle the 
most simple model describing with acceptable precision the evolution of East- 
erner society. With the idea of writing a model that could be generalized to 
a more complex system, we first divide the elements into two categories: the 
resource quantity Ri with i — 1, 2 • • • k and the inhabitants numbers (species) 
Ni with % = 1, 2 • • • m. With the concept of resources we are meaning resources 
in a very large sense, it could be oil, trees, food and so on. The several kind 
of resources are described by the index i. In similar way, with the concept 
of inhabitants, we are meaning different species of animals or internal sub- 
division of human being in country or town or even tribes. Leaving the idea 
of a constant quantity of resources that leads to the logistic equation for the 
number on inhabitants [7], the aim of this paper is to include in the dynamical 
description of the time evolution of the system the resources too which can 
not be considered as constant. A generalization of the logistic equation to an 
arbitrary number of homogeneous species interacting among the individuals 
(with non constant resources) can be written as 

d T N- 1 m 

—Ni = ri Ni 1 - Y XaNjNi. (1) 

dt 1 [ N ci (R l ---R k )\ j= ^j J 3 K ' 



Where is the usual growing rate for species i. In the denominator it appears 
the carrying capacity of the system with respect to the number of inhabitants 
N C i (Ri ■ ■ ■ Rk). Beside the dependence on the resources Ri we could have also 
a dependence on an other species that would be then a "resource" for some 
other species. This fact is expressed even by the quantities Xij that in general 
are not symmetric expression (xij Xji) since that the prey is a resource for 
the predator and not the inverse. Similarly, for the resources we have: 

m 

-^(XijNjRi. (2) 
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It is clear that set of equations ([2]) could be formally included into Eq. ([T]) 
redefining the quantities a^. Nevertheless, we shall keep this distinction for 
the sake of clarity especially referred to resources, such as trees, oil or oxygen, 
where the carrying capacity is not determined by other species and can be 
considered a constant (R C i). Also the meaning of the parameters r[ is more 
or less the same of the analogous parameters ?v It is suitable to define r[ as 
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" renewability ratios", since describe the capacity of the resources to renew 
itself and clearly are depending on the kind of resource. For example, the 
renewability of the oil is clearly zero since the period of time to get oil from 
a natural process is of the order of geological time-scale processes. In general 
all parameters of Eqs. (JT]) and © are time dependent, including stochasticity. 
We can assume reasonably slowly time-varying for the ancient societies so 
that we can consider it as constant [8], particularly the a's. Anyway the set of 
cty is worthy of a more detailed discussion. We can call this set of parameters 
technological parameters in the sense that they carry the information about the 
capacity to exploit the resources of the habitat. We shall see in the next section 
that the technological parameters combined with the renewability ratios will 
be the key point to decide wether, or not, a society is destined to collapse. 



2 Easter island collapse model 



The particular history of Easter Island society presents several advantages for 
modelling its evolution [9J. In fact it can be with very good approximation 
considered a closed system. The peculiar style of life and culture allows us to 
consider a basic model where trees are essentially the only kind of resources. 
Many of activities of the ancient inhabitants involved the trees, from building 
and transport the enormous Moai, to build boats for fishing, etc. In fact the 
cold water was not adapt to the fish life and the impervious shape of the coast 
made difficult fishing. Finally from historical reports it can be inferred that 
the inhabitants did not change the way to exploit their main resource, even 
very near to exhaust it so that we can consider the technological parameter 
as constant. Considering Eqs. (pQ) and (T5]) for one inhabitant species and one 
kind of resource, we obtain: 



d 
~dt 



N = rN 



1 - 



N 



N C (R) 



—R = r'R 
dt 



R c 



aNR, 



(3) 
(4) 



where we introduce the notation: «u = a. The unknown function N c (R) 
has to satisfy few properties. For a quantity of primary unlimited resource, 
R — > oo, even iV c (R) — > oo that means that the population can grow unlimited 
too. In the opposite case R — > clearly also the population must vanish, 
N c (R) — > 0, and finally when the resource is constant we are back to the 
ordinary equation (logistic) so that N c (R) = const. It is clear that the choice 
of this relation is quite arbitrary but following the simplicity criteria we can 
select N c (R) = (3R, where j3 is a positive parameter. This choice formalizes 
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the intuition that the maximum number of individuals tolerated by a niche is 
proportional to the quantity of resources. We note that in (j3J) the interacting 
term depend on the variable R. Namely, for R = no variation of resource 
exist {4rR = 0) corresponding to a biological criterion and different from this 
one of reference [TP]. A more sophisticate model should include also fishing as 
resource and consider for N c (R) an expression such as N c (R) = (3R + Rf 
where Rf is the fishing carrying capacity This resource was limited near the 
coast and could not be fully exploited without boats, so that, we are going to 
neglect this resource. We can rewrite Eqs. ([3]) and (T4|) as: 

N] 

'-ml' (5) 

where oce = — f . (6) 
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The dimensionless parameter is the ratio between the technological pa- 
rameter a, representing the capability of to exploit the resources, and r' the 
renewability parameter representing the capability of the resources to regen- 
erate. We will call a E deforestation parameter since it gives a measure of the 
rapidity with which the resources are going to exhaust and then a measure of 
the reversibility or irreversibility of the collapse. Using the historical data we 
can have an estimation of the parameters. At the origin (t = 0) we can assume 
that the trees were covering the entire island surfaces of 160 km 2 . When the 
first humans arrived to the island, around the 400 A.D., their number was of 
the order of few dozens of individuals and it grew until to reach the maximum 
N M ~ 10000 around the 1300 A.D. 

Finding the equilibrium points of Eqs. ((5]) and fl6j we obtain (see section III, 
for stability): 

iVo = 0, Ro = R c (unstable-saddle-point) (7) 

e l + otEfiRc e I + OLePRc ^ ^ ^ 



While the point (N , R ) of Eq. ([7]) represents the trivial fact that in absence 
of human being the number of trees is constant (carrying capacity). Eq. 
describes the fact that, due to the interaction humans-environment, the more 
interesting equilibrium point (N e ,R e ) does not coincide with (N C ,R C ) with 
N c = (3R C since a ^ 0. To study the stability of the point (N e ,R e ), we have 
linearize the system of Eqs. ([5]) and (J6j around the equilibrium point. In fact, 
in the next section we shall show that it is a stable equilibrium point. 
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3 Equilibrium points and stability 



Let us first cast Eqs. (jSJ) and (jSD in term of dimensionless quantities; setting 
u(t) = N/N c , g(t) = R/R c and r = rt, we have: 



-v = v 



d , v 

dr 

d r' 
— g — rg [1 — g — ocu] , a = oeN c , r = — 
dr r 



(9) 
(10) 



For sake of clarity we rewrite also the equilibrium point: 

1 1 

= —— , Qe = — — 
1 + a 1 + a 



with obvious meaning of the symbols. Perturbing the equilibrium point v(t) = 
v e [1 + T}(t)] and g(r) = g e [1 + s(r)] with r/(r) and e(r) infinitesimal func- 
tions, after straightforward algebra we obtain Eq. (jSJ): 



d 

d7 f 



v 



-rj + e 



1 + a 



[arj + e] 



:i2) 



(13) 



The eigenvalues of the system are: 



Al,2 



-(1 + a + f) ± \J(l + a + f) 2 - 4f(l + a) 
2(1 + a) 



(14) 



Restricting ourself to the case of positive values of a, Eq. (1141) shows that both 
the eigenvalues always have a real negative part, so that the equilibrium point 
(v e , g e ) is a stable equilibrium point. More in detail we have that for 



1 r _ 1 _ ( VT ~ 1 

r < -, a > and for r > -, a < - — 1= — 
~ 4' - 4' - 2-v/f -1 



the eigenvalues are real and negative so that the equilibrium point is reached 
in exponential damped way, otherwise the eigenvalues acquire an imaginary 
part and the system reach the equilibrium point via exponential damped os- 
cillations. 
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Fig. 1. N(t)/N c versus time. In this numerical example N(0)/N c ~ 0.8 and 
N Max /N c ~ 3 

4 Collapse condition 

Even if, mathematically speaking, the stable point (N e , R e ) is an acceptable 
result we have to take in account the biological constraints that allow to a 
specie to survive. A reasonable number of individuals is required for viability 
of a given species [TT][T2] . This is so because genetic diversity, social structures, 
encounters, etc., need a minimum numbers of individuals since under this 
critical numbers the species is not viable and collapse. It is worthy to stress 
that while the trees can reach an equilibrium point without the humans, Eq. 
(J7|), the opposite does not hold, as stated by Eq. (jHJ). 

Calling the minimum number of humans N m i n we can find a upper bound for 
the parameter «e so that a civilization can survive. Imposing the condition 
that at the equilibrium point N e > N min we obtain: 

ob > — — — (collapse condition). (15) 

-'' min P-ft'C 

As further simplification of inequality ffl~5]) . we assume that N min <C (3R C and 
we find that aEN min > 1 or aN min > r' . It is a natural condition since it tells 
that collapse exists when the production rate r' is minor than the deforestation 
rate aN min . 

More in general it can be showed, numerically, that considering the standard 
case with the starting population number N(0) < N c = (3R C , we can have a 
solution that can exceed the value /3R C (depending on initial conditions N(0) 
and T^| t= o). It is worthy to stress that in the case of ordinary logistic map 
the population number never can exceed this limit value|7j. In the example 
of Fig. (U the maximum of N is reached at a value that is almost three times 
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f3R c . Then, according to the region of the parameters that we are considering, 
the paths to the final equilibrium is exponentially fast, reaching eventually the 
point N min and collapsing. 



5 Deforestation rate estimation 



As we saw in the previous section, the collapse condition (15) gives a sufficient 
condition on the deforestation rate per individual a . On other hand, the last 
period of tree extinction was governed essentially by the deforestation rate. In 
this way, we have the rate of tree extinction 

ldR 

-- ~ -alt. (16) 



As discussed, the path to the equilibrium point is exponentially fast, so that a 
rough estimation of the left side of Eq. ( lT6"j) is the time scale of the deforesta- 
tion, t>, while the right side can be taken at the end of the collapse process 
(the equilibrium point): 

— ~ aN F , (17) 



being Np the final number of individuals. It can be deduced [13] that the range 
of is Tp ~ 100 yrs. to Tp ~ 300 yrs. and Np ~ 3000, the rate of deforestation 
(per individual) could be estimated as: 

a ~ ( yrs. individual) -1 (18) 

TpNp 



giving a range 

1.1 x 10~ 6 < a < 3.3 x 10 -6 ( yrs. individual)' 1 . (19) 



This estimation has validity in the case of exponential decay which is the our 
case, as it has been showed by the analysis performed in Sec. [3j Assuming 
the number of trees as proportional to the area A we can now estimate the 
rate-deforestation-area. The island has a surface at order of A ~ 160 km 2 and 
initially it can be supposed that was covered of trees so that we can estimate: 

dA An , 9 , . , 

0.5 < < 1.6 (km 2 /yrs). (20) 

dt Tp 
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As comparison we can consider that in the last 500 years the deforestation 
of amazonian forest rate is 15 (km 2 /yrs). Considering that in 500 years the 
deforestation technology became more and more efficient, especially in the 
last century, we can consider it as an upper limit, giving us an idea of the 
technological change. In the human history there are several examples of over 
exploiting the natural resources even if not so known as Easter island. In 
particular, the Copan Maya history[T3] has certain similarity with respect to 
the technology level and the over exploiting of the natural resources. In short, 
this ancient civilization reached almost 20.000 individuals and declined to 
5000 individuals in the 9th century. Using the estimation of the technological 
parameter a obtained for eastern island civilization, we get a collapse time 
from Eq. (fl7|) that is 60 < r < 180 years. The collapse time based on historical 
reports is r ~ 100, showing that the adopted model is consistent with the 
available data. The estimation of the parameter is consistent with the idea 
that similar civilizations, in technological meaning, have similar capacity to 
exploit the natural resources. 



6 Concluding remarks 

A mathematical model considering the interaction among carrying capacity 
and population in an isolated system has been considered. The model takes 
in account the fact that a population can over exploit the carrying capac- 
ity without saturate, a fact of relevant importance on the path leading to a 
collapse of a society. Its application to the collapse of the Easter Island civi- 
lization has been presented. An estimation of the technological parameter a is 
obtained and applied to an other ancient civilization, the Copan-Maya, with 
a reasonably precise expectation about their collapse time. All confirming the 
consistency of the adopted model. On the other hand, its relative reasonable 
prediction suggests a possible extension to more complex system. The effort 
to mathematical modelling of ancient civilizations [6f8fT5] could be important 
considering the actual human growing and resources exploitations. An ade- 
quate equilibrium between competition, demand and exploitation is the key 
of surviving. 
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